home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zpptrs.z / zpptrs
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. ZZZZPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          ZZZZPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZPPTRS - solve a system of linear equations A*X = B with a Hermitian
  10.      positive definite matrix A in packed storage using the Cholesky
  11.      factorization A = U**H*U or A = L*L**H computed by ZPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE ZPPTRS( UPLO, N, NRHS, AP, B, LDB, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDB, N, NRHS
  19.  
  20.          COMPLEX*16     AP( * ), B( LDB, * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      ZPPTRS solves a system of linear equations A*X = B with a Hermitian
  24.      positive definite matrix A in packed storage using the Cholesky
  25.      factorization A = U**H*U or A = L*L**H computed by ZPPTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangle of A is stored;
  31.              = 'L':  Lower triangle of A is stored.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      NRHS    (input) INTEGER
  37.              The number of right hand sides, i.e., the number of columns of
  38.              the matrix B.  NRHS >= 0.
  39.  
  40.      AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
  41.              The triangular factor U or L from the Cholesky factorization A =
  42.              U**H*U or A = L*L**H, packed columnwise in a linear array.  The
  43.              j-th column of U or L is stored in the array AP as follows:  if
  44.              UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j; if UPLO =
  45.              'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
  46.  
  47.      B       (input/output) COMPLEX*16 array, dimension (LDB,NRHS)
  48.              On entry, the right hand side matrix B.  On exit, the solution
  49.              matrix X.
  50.  
  51.      LDB     (input) INTEGER
  52.              The leading dimension of the array B.  LDB >= max(1,N).
  53.  
  54.      INFO    (output) INTEGER
  55.              = 0:  successful exit
  56.              < 0:  if INFO = -i, the i-th argument had an illegal value
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.